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ABSTRACT: We study the differential branching ratio, branching ratio and the forward-backward 
asymmetry for the exclusive B — ► ir£ + £~md B — ► p£ + £~ decays in the general two Higgs doublet 
model including the neutral Higgs boson effects. We analyze the dependencies of these quantities 
on the neutral Higgs boson contributions and the other model parameters. We observe that two 
Higgs doublet model with the neutral Higgs boson exchanges gives quite sizable contributions to 
these observables for both channels we consider. Since the neutral Higgs boson exchanges are the 
only source of the forward-backward asymmetry for B — ► 7rT + r~decay, which is at the order of 
magnitude 1 — 10%, measurement of this observable is promising to determine the neutral Higgs 
boson effects. 
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1. Introduction 

The rare decays of B-mesons, induced by the flavor-changing neutral currents (FCNC), have always 
been a good candidate for testing the Standard Model (SM) at the loop level and looking for new 
physics beyond it. They can also be used to determine the fundamental parameters of the SM, like 
the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, the leptonic decay constants etc.. 
Among the rare B -decays, exclusive processes induced by b — > s{d)£ + £" transitions have received 
a special attention since the SM predicts relatively larger branching ratios for these decays. 

From the experimental side, there is an impressive effort for searching B -meson decays, espe- 
cially at the B-factories, like at Belle [jl|] and BaBar [0], and with the increased statistical power 
of these experiments, in the near future the rare B-meson decays will be measured very precisely. 
From the theoretical side, the decays B — > X s ^£ + £~ provide important probes of the effective 
Hamiltonian which governs the FCNC transitions b — > s(d)£ + £~ at quark level [[J. For b — > s£ + £~ 
transition, the matrix element contains the terms that receive contributions from tt, cc and uu loops, 
which are proportional to the combination of £ t = VtbV t * s , £ c = V c bV* s and £ u = V u bV* s , respec- 
tively. Smallness of ^ u in comparison with £ c and together with the unitarity of the CKM matrix 
elements, bring about the consequence that matrix element for the b — > s£ + £~ decay involves 
only one independent CKM factor ^, so that the CP violation in this channel is suppressed in 
the SM [§ |]. However, for b d£ + £~ decay, all the CKM factors rj t = V tb V* d , r/ c = V cb V* d 
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and rj u = V u t,V* d are at the same order in the SM and this induces a considerable CP violating 
asymmetry in the partial rates |§, 0]. 

In this paper we investigate the exclusive B — > ir£ + £~m& B — ► p£ + £~ decays, which are 
induced by the b — ► d£ + £~ decay at the quark level, in the framework of the general two Higgs 
doublet model (2HDM) (model III). These decays have been studied in the literature, both in the 
SM and in the 2HDM. CP violating effects in inclusive b — ► d£ + £~ and exclusive B — ► 7r^ + £~and 
B — > p£ + i?~channels were studied within the framework of the SM in refs. [§]-[§]. 2HDM 
contributions to these exclusive decays have been investigated in |ffi,|IOp. In earlier works about the 
exclusive B — > p£ + £~and B — ► 7r£ + ^decays, contributions from exchanging the neutral Higgs 
bosons (NHB) were neglected because of the smallness of m^/mw {£ = e, //). However, in the 
models with two Higgs doublets, such as MSSM, 2HDM, etc., the situation is different, especially 
in case of £ = r with large tan (3 or large C/^ TT , where tan/3 = V2/V1, the ratio of the vacuum 
expectation values of the two Higgs doublets and ^ rr is the Yukawa coupling to which NHB 
contributions are proportional in model III. Indeed, there are a number of works in the literature, 
which show that the contributions from exchanging NHB can compete with those from exchanging 
7 and Z when tan/3 and/or TT are large enough [11 ]-[|l6|]. 

The aim of this work is to calculate the branching ratio (BR) and the forward backward asym- 
metry (Afb) of the exclusive B — > -K£ + £~mA B — ► p£ + ^ _ decays in the general 2HDM, including 
NHB effects without neglecting the lepton mass. The 2HDM is one of the simplest extensions of 
the SM, which is obtained by the addition of a second complex Higgs doublet. In general, the 
2HDM possesses tree-level FCNC that can be avoided by imposing an ad hoc discrete symmetry 
[17]. As a result, there appear two different choices, namely model I and II, depending on whether 
up-type and down-type quarks couple to the same or two different Higgs doublets, respectively. 
Model II has been more attractive since its Higgs sector is the same as the Higgs sector in the su- 



persymmetric models. In a more general 2HDM, namely model III [|18|, |19p, no discrete symmetry 
is imposed and there appear FCNC naturally at the tree level. We note that in model III, FCNC 
receiving contributions from the first two generations are highly suppressed, which is confirmed by 
the low energy experiments. As for those involving the third generation, it is possible to impose 
some restrictions on them with the existing experimental results. 

The paper is organized as follows: In section ||[ after we present the theoretical framework 
of the general 2HDM and the leading order QCD corrected effective Hamiltonian for the process 
b — ► d£ + £~ , we calculate the differential BR and the Afb of the exclusive B — ► p£ + £~&n& 
B — ► n£ + £r decays. The ||. section is devoted to the numerical analysis and the discussions. 
Finally, in the Appendices, we give the explicit forms of the operators appearing in the effective 
Hamiltonian and the corresponding Wilson coefficients. 



2. The exclusive B — > n£ + £ and B — > p£ + £ decays in the framework of the general 
2HDM 

In this section, we first present the theoretical framework of the general 2HDM, and then calculate 
some physical observables related to the exclusive B — * ir£ + £^a.nd B — > p£ + £~ decays in this 
model. 
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2.1 The theoretical framework 



We would like to present the main essential points of the general 2HDM, namely model III. In this 
model, both Higgs doublets can couple to up- and down-types quarks. We can choose two scalar 
doublets 4>\ and 4>2 in the following form 



01 = 7I 
with the vacuum expectation values, 




+ 



J_ / V2H+ 
^2 \ Hi + iH 2 



(2.1) 



<4>i>=—=\ ; < 02 >= 



(2.2) 



wo 

so that the first doublet <j)\ is the same as the one in the SM, while the second doublet contains all 
the new particles. Further, we take H\ and H 2 as the mass eigenstates h° and A , respectively. 
The general Yukawa Lagrangian can be written as 



-) - VijQiKpiUjR + rifjQiL(t>iDj R + £ y ^QiLfoUjR + ZfjQiLfaDjR + h.c. 



(2.3) 



where i,j are family indices of quarks , L and R denote chiral projections L(R) = 1/2(1 ^75), (f>m 
for m = 1,2, are the two scalar doublets, Qn are quark doublets, UjR, DjR are the corresponding 
quark singlets, rjfj D and £^- D are the matrices of the Yukawa couplings. After the rotation that 
diagonalizes the quark mass eigenstates, the part of the Lagrangian that is responsible for the FCNC 
at the tree level looks like 



C y ,fc = ~H ] U [Vckm $ f Vckm L] V , (2.4) 

where U(V) represents the mass eigenstates of up (down) type quarks. In this work, we adopt the 
following redefinition of the Yukawa couplings: 



JJ.D 
?7V 



4Cf tU,p 
V2 ^ N ' ir 



(2.5) 



After this brief summary about the general 2HDM, now we would like to present briefly the 
main steps in calculating the matrix elements for the inclusive b — ► d£ + £~ decay. For this, the 
effective Hamiltonian method provides a powerful framework. In this approach, the first step is to 
calculate the full theory including the NHB effects. We use the on-shell renormalization scheme 
to overcome the logarithmic divergences that appear during the calculations of NHB contributions. 



(For the details of this calculations, see ref.[]l The next step is to match the full theory with the 
effective theory, which is obtained by integrating out the heavy degrees of freedom, i.e., t quark, 
W^, H^, h°, H° and ^4° in our case, at high scale /x = mw, and then calculate the Wilson 
coefficients at the lower scale n ~ 0(rrib) using the renormalization group equations. Following 
these steps above, one can obtain the effective Hamiltonian governing the b — ► d£ + £~ transitions, 
in the 2HDM in terms of a set of operators 



H 



eff 



4G F 

~7f 



v tb v* d 




10 



i=i 



K{Ci(^)[ou^) - o 1 (p)] + c 2 {n)[o u M - o 2 (m)]} 



(2.6) 
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where 

using the unitarity of the CKM matrix i.e. V t bV t * d + V ub V* d = -V cb V* d . Here, Oi and 2 are 
f/ie current-current operators, Oz,...,0§ are usually named as the QCD penguin operators, Oj and 
Og are ?/ze magnetic penguin operators and Og and O10 are semileptonic electroweak penguin 
operators. The additional Qj (1 = 1, 10) are due to the NHB exchange diagrams. Of-O^ are 
the new operators for the b —> d decay which are absent in the b —> s decay. Cj(ju) and Cq 1 (p) are 
Wilson coefficients renormalized at the scale p. All these operators and the Wilson coefficients, 
together with their initial values calculated at \i = mw in the SM and also additional coefficients 
coming from the new Higgs scalars are presented in Appendix |A|. 

Neglecting the mass of the d quark, the effective short distance Hamiltonian for the b — ► d£ + £~ 
decay leads to the QCD corrected matrix element: 

M = ^^V tb V t * d ^C e 9 ff (m b ) J 7m (1 - 75)^7^ + C 10 (m b ) J 7/i (l - l5 )bl^ l5 £ 

- 2C e 7 ff (m b ) ^dia^il + lh )bhn + C Ql (m b ) d(l + 75)6^ + Cg 2 (m b ) d(l + 7s)W 



(2. 



where g is the momentum transfer. 

2.2 The exclusive £ n£ + £~ decay in the 2HDM 

In this section we calculate the Si? and the ^4fb of the B — > 7r^ + ^~decay. In order to find these 
physically measurable quantities at hadronic level, we first need to calculate the matrix elements 
< 7r(p*)|d 7/x (l - 7 5 )&|5(p jB ) >, < 7^)^07^(1 + 75)&|^(pb) > and < ir(p w )\d(l + 
75)6|B(pb) >. The first two of these matrix elements can be written in terms of the form factors 
in the following way 



(2.9) 



< Tr(p 7V )\dia^q u (l+ -y 5 )b\B(p B ) > = \{pb +Pir)f l q 2 - q^{m 2 B - ml)]f v (q 2 ) , (2.10) 

where p B and p^ denote the four momentum vectors of B and 7r-mesons, respectively. To find 
< ir(p n )\d(l + ^5)b\B(p B ) >, we multiply both sides of Eq. with and then use the 

equation of motion. Neglecting the mass of the <i-quark, we get 

1 



< 7r( Pn )\d(l + 75 )b\B( PB ) >= _[/+(^)(™b - mi) + f-{q 2 )q Z ]- 

m b 



(2-11) 



Using Eqs. ( |2.9r|2.11| ), we find the amplitude governing the B — ► 7r£ + £ decay : 

sB^TT _ GpOL 



2\[2ix 



V tb V*A [2A lP £ + B iq »]£^£ + [2G lP £ + D iq »]£ W £ + E x ££ + F x £ lh 



(2.12) 
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where 



Ax = C e 9 ff f + - 2m B C e 7 ff f v , 
B 1 =C e 9 ff (f+ + f-)+2C e 7 ff ^f v ( 

G\ = Cio/ + , 

Di = c 10 (/ + + r), 

E 1 =C Ql j-[(m 2 B -ml)f+ + f- q % 
F 1 = C Q2 ±-l(m 2 B -ml)f + + f-q 2 }. 
Here f + , f~ and f v are the relevant form factors. 



2 2 2\ 

m B -m v -q ) 



(2.13) 



The matrix element in Eq. ( |2.12| ) leads to the following double differential decay rate: 
<fiY __ G^ 2 , W |V^ v\m%)*{l-W)\A x \* + atf\Etf + \F x \*) 



ds dz 2 n 7r 5 



+ (m|A 7r (l - zV) + 16 r w mj) \Gi\ 2 + 4 s mj |£>i| 2 



+ 4 mf (1 - rv - a) i2e[GiDJ] + 2 u m £ z i2e[Aif;j] 
+ 2 m £ ((1 - r w - a) i?e[GiF*] + sRe[D x Ft}) \ . 



(2.14) 



Here s = q 2 /m B , r w = m 2 /m 2 B , v = J 1 — i = mi/mB, K = r 2 + (s — l) 2 — 2r,r(s + 1) 
and 2 = cos 9, where 9 is the angle between the three-momentum of the £~ lepton and that of the 
B-meson in the center of mass frame of the dileptons £ + £~ . We note that our expression for double 



differential decay rate in Eq. (2.14) coincides with the one in [2C]. 



Integrating the expression in Eq. ( |2.14| ) over the angle variable, we obtain the differential 
decay rate as follows 



dr 

ds 



where 



A, 



Q2 a 2 

^5 I V tb V t * d 1 2 m 3 B ^/K v A w , 



4m 2 



(2.15) 



m B X n (3 - ^Xl^l 2 + |Gi| 2 ) + -5-^(12 r w s + X^Gx 



3s 



+ 4mjs |L»i| 2 + s(v 2 |£i| 2 + |Fi| 2 ) + 4 mj(l - 
+ 2 m/((l -r n - s)Re[G 1 F*] + sRe[D x Ff]) 



s)Re[G 1 Dl\ 



(2.16) 



The Afb is another observable that may provide valuable information at hadronic level. We 
write its definition as given by 



A F b(s) 



(2.17) 



where F is the total decay rate. 
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For the B — > ir£ + £ decay, Ap B is calculated to be 

A FB = - j ds {tv 2 KRe{Ai ED) I J ds v^/)^A n . (2.18) 



As seen from Eq.(2.18), in the B — > ir£ + £ decay, the only source for the Ap B is the NHB effects 



[12, EQ]. Since A FB does not exist in the SM and 2HDM without NHB effects, it may be a good 



candidate for testing the existence and the importance of the NHB contributions . 
2.3 The exclusive B -> pl+f decay in the 2HDM 

In this section we proceed to calculate the BR and the Ap B of the B — > /9£ + ^ _ decay. We follow 
the same strategy as in the B — ► 7r^ + ^~ decay. In order to calculate the matrix element governing 
the — > p£ + ^ _ decay, we need the following matrix elements: 

2V(v 2 ) 

< p( Pp ,e)\d^(l - j 5 )b\B(p B ) > = -e puXa e* u p x p B ^-L- - ie*(m B + m^A^q 2 ) 

F m B + m p p 



+ i(PB +Pp)^*Q) ^ Q 1, + iq^eq)^[A 3 {q 



p r /W„2> 

A)(<Z 2 )], (2-19) 



4e^ A(7 £ > >yT 1 (< Z 2 ) + 2i[e* (m| - m 2 ,) 
(PB+pA*(eW2(? 2 ) + 2i(e*g) 

/ «2 \ 



9m - (PS + £>A 



m 



r 3 (^), (2.20) 



< p(p p ,e)\d(l + j 5 )b\B{p B ) >= —2im p (eq)A 



(2.21) 



where p p and e* denote the four momentum and polarization vectors of the p meson, respectively. 



In order to calculate the matrix element in Eq. ( E.21[ ), we multiply both sides of Eq.( E.19[ ) with q p 
and use the equation of motion. 



From Eqs. (2.19-2.21), we get the following expression for the matrix element of the B 

p£ + £~decay: 



M 



B^p 



Gfcz 



„ ^WjVPV^ + iBe; - iC(p B +P P Me*q) - iD(e*q)q p ] 
+ h^l[1E^ uXfT e* v p x p p B + iFel - iG{e*q)(p B + p p ) - iH{e*q)q p ] + i££Q{e*q) 
+ i£j 5 £N(e*q) I (2.22) 
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where 



A = C?'—?—+4?±C?'T 1 , 
ms + m p q z 



B = (m B + m p 



\C e / f A 1 + ^(m B - m p )C e 7 Si T^ , 

C = C^^ 2 — + 4^C?" U + -^— 2 T^\ , 
m B +m p q z \ m B — m z J 

D = 2C e 9 ff ^(A 3 - A ) - 4C?^T 3 , 

E = Cio - , (2.23) 

F = C 10 (m B + m p )Ai, 
G = C W - M 



m B + m p 

2 ' ' A * 



H = 2C 10 ^(A 3 -A ), 



T 

Q = 2Cq 1 — -A , 
m b 

N = 2C Q2 ^A . 
m b 

Here Aq, A\, A%, A 3 , V, T\, T 2 and T 3 are the relevant form factors. 

This matrix element leads to the following double differential decay rate 

d 2 T a 2 G 2 



(]s ( ,_ , >r , m , )jT 5 \ v tbV t d\ 2 y%v I 4 s \ p {2 + v 2 (z 2 - l))\A\ 2 + 4v 2 s m%X p {l + z 2 )\E\ : 



+ 16 ml sv z,/\~J Re\ BE*} + RelAF*}) +- 



[A p (l - zV) + 2r p s(5 - 2v 2 )]\B\' 



+ m%X 2 p (l - z 2 v 2 )\C\ 2 + [A p (l - z 2 v 2 ) - 2r p s{\ - Av 2 )]\F\ 2 
+ m%X p [{-l + r p ) 2 (l - v 2 )z 2 + (-1 + z 2 )(st 2 - 8(1 + r p )t 2 - X p )]\G\ 2 
+ 2m 2 B X p W p (l - z 2 v 2 )Re[BC*} - 2m 2 B X p [W p {\ - z 2 v 2 ) - At 2 }Re[FG*} 
+ m 2 B X p (4sm e (m e \H\ 2 + Re[HN*]) + s(\N\ 2 + v 2 \Q\ 2 ) - 4t(Re[F{2tH* + N*/m B )} 



+ 4(1 - r p )m e Re[G(2m e H* + N*)]j + Um B vz z Re[(W p B + m B {W z p - Ar p s)C)Q*] 

(2.24) 

where r p = m 2 /m 2 B , W p = — 1 + r p + s and X p = r 2 + (s — l) 2 — 2r p (s + 1). 
The differential decay rate for B — > p£ + £~ decay is given by 

dT a 2 G 2 F m B ,„ T/ * |2 /r- A , OCk 

rfj = 2 i 27r 5 V*P w A p (2-25) 

where 

A p = ^A p m|s((3 - w 2 )|^| 2 + 2v 2 |£| 2 ) - -A p m|m^e[(F - m|(l - r p )G - m 2 B sH)N*\ 
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+ — \ p m B 



sv 2 \Q\ 2 + \x p m 2 B (3 - v 2 )\C\ 2 + s\N\ 2 + m 2 B s 2 (l - v 2 )\H\ 2 



+ oK 3 - v > W P - 3 s ( x - v 2 )]Re[F G*] — 2 s (1 — v )Re[F H*] 



+ 2m|s(l - r p )(l - u 2 )i?e[G #*] + -(3 - v 2 )W p Re[B C*} 

3 



1 2 

3r n 



(A p + 12r p s)(3 - v 2 )|£| 2 + A p m|j[A p (3 - v 2 ) - 3s(s - 2r p - 2)(1 - i> 2 )]|G|< 



+ (A p (3-v 2 ) + 24r p st; 2 )|F| 2 



(2.26) 



We also give the Afb of the 5 — ► decay 



,4 Fi j = J ds2m B \ p v 2 (4m B s(Re[B E*]+Re[AF*]) 



+ —[W p Re[B Q*] + m 2 B \ p Re[C Q*}} 



dSyJXp V Ap. 



(2.27) 



3. Numerical results and discussion 



In this section we present the numerical analysis of the exclusive B — > irl + £~a.nd B — > p£ + ^~decays 
in the general 2HDM. We give our results for the £ = t case in order to express our motivation 
to calculate the BR and the Afb of these decays without neglecting the lepton mass. The input 
parameters we used in our numerical analysis are as follows: 

m B = 5.28 GeV , m b = 4.8 GeF , m c = lAGeV , m T = 1.78 GeV, = 0.14 GeF, 

m p = 0.77 GeF , m H o = 150 GeF , m h o = 100 GeF , m A o = 100 GeF , m H ± = 400 GeV , 



\V tb V t * d \ = 0.011 , a" 1 = 129 , G F = 1.17 x 10~ 5 GeF~ 2 , r B = 1.54 x 10" 12 s . 



(3.1) 



Using the Wolfenstein parametrization of the CKM matrix, \ u in Eq.([2.7[) can be written as: 

p(l — p) — r/ 2 — if] 



A t( 



(1 - p) 2 + r] 2 



+ 0(A 2 



(3.2) 



Furthermore, we have used the relation 

\VtbV t *f _ 



A 2 [(l- pf + ^j+OiX 4 



(3.3) 



where A = sm9c — 0.221 and we take the Wolfenstein parameters as p = —0.07 and rj = 0.34 
throughout the calculations. 

We note that the Wilson coefficient Cg receives also long distance (LD) contributions due 
to cc intermediate states. (See Appendix ||] for the details of the LD contributions). 
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There are five possible resonances in the cc system that 
can contribute to the decays under consideration and to calcu- 
late their contributions, we need to divide the integration re- 
gions for s into two parts: the region 0.4546 < s < ({m^ 2 — 
0.02) /tub) 2 is common for both decays and we have ((m^ + 





7(0) 


Tf 


n f 




0.29 


6.71 


2.35 


f B->« 


-0.26 


6.553 


2.30 


£B — >7T 
J V 


-0.05 


6.68 


2.31 



0.02)/m B ) 2 < s < 0.9476 and ((m^ + 0.02)/m B ) 2 < s < 
0.7296 as second integration parts for the B — ► 7rr + r~and 
B — » /rr + T~decays, respectively. Here, m^ 2 is the mass of quark model 
the second resonance. 

The masses of the charged and the neutral Higgs 
bosons, m H ±, m A o, m h o, m H o and the Yukawa cou- 
plings ) remain as free parameters of the model. 
For the mass of the charged Higgs, the lower limit 
m H+ > 200 GeV and m H+ > 250 GeV have been 



Table 1: B — *■ it transition form 
factors in the light-cone constituent 







F(0) 


CLF 


b F 






0.26 ± 0.04 


0.29 


-0.415 


At 


-*p 


0.22 ±0.03 


0.93 


-0.092 


yB- 




0.34 ±0.05 


1.37 


0.315 


T B- 


^p 


0.15 ±0.02 


1.41 


0.361 


Tt 


->p 


0.15 ±0.02 


0.28 


-0.500 


T 3 B ~ 


+p 


0.10 ±0.02 


1.06 


-0.076 



given in [ |22| ] and [|23|], respectively. However, it is 
also pointed out that adding different theoretical errors 
leads to mjj ± > 370 GeV and also these bounds are 
quite sensitive to the details of the calculations. In our 
work, we choose m# ± = 400 GeV . For the masses of 
the neutral Higgs bosons, the lower limits are given as 

m H a > 115 GeV, m h o > 89.9 GeV m A o > 90.1 GeV in p4| ] and the values we choose are given 
in Eq.(3.1). 



Table 2: B — ► p transition form factors in a 
three-parameter fit. 



For Yukawa couplings , we use the restrictions coming from CLEO data [[25]], 
BR{B X s 7) = (3.15 ± 0.35 ± 0.32)10" 4 , 



(3.4) 



B 



o 



B mixing [|26|], p parameter 



J, and neutron electric-dipole moment []27[], that yields 



^Nib ~ *-* an( ^ ^Nij ~ 0, where the indices i, j denote d and s quarks, and £jy tc << ^ tt - 
Therefore, we take into account only the Yukawa couplings of b and t quarks, ^ tt , ^Nbb an ^ a l so 



£jy TT . There is also a restriction on the Wilson coefficient Cy from the BR of B 



Eq.(gJ) as follows U 



0.257 < \C? ff \ < 0.439. 



X s j in 



(3.5) 



In the following subsections, we calculate the dependencies of the differential branching ratio 
(dBR/ds) and Afb(s) of the above decays on the invariant dilepton mass s, and also dependencies 
of the BR and Afb on the Model III parameters. The results are presented by a series of graphs, 

I 8 



which are plotted for I 
or r th > 1. 



r and for two different cases of the ratio |r t b| 



CD 



where \r t b\ < 1 



3.1 Numerical results of the exclusive B — > ir£ + £ decay 

In our numerical calculations for B — » vr^+^decay, we use the results of the light-cone constituent 
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quark model [28, 29], which can be found from the following expression 

/tf) = m 

{l-q 2 /T f ) n f' 
where the parameters /(0), Tf and nj are listed in table [I]. 



(3-6) 




0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 

s 



Figure 1: dBRj ds for B — ► 7TT + T~as a function of s for bb = 40 mt and £^ rT = 10 m T , in case that 
the ratio \r t b\ < 1. Here the region between the solid curves represents the dBR/ds in Model III without 
the NHB effects, while the one between the dashed curves is for the dBR/ds with NHB contributions. The 
SM prediction is represented by the small dashed curve 
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Figure 2: The same as Fig.dU), but for r t b > 1 with bb = 0.1 mj and £® TT = m 



In Fig. ([T]), we plot the dependence of the dBR/ds on the invariant dilepton mass s, for 
\r tb \ < 1 and C e 7 Sf > cases, by taking into account the long distance effects. Here the region 
between the solid (red) curves represents the dBR/ds in Model III without the NHB effects, while 
the one between the dashed (green) curves is for the dBR/ds with NHB contributions. In both 
cases, we observe an enhancement compared to the SM prediction, which is represented by the 
small dashed (blue) curve. This enhancement reaches up to 70 % and 50 % as compared with 
the SM and the Model III prediction without NHB effects, respectively. In Fig. (g), the same 
comparison is made for r#, > 1 case and we see that the contributions coming from the NHB effects 
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40 50 

^N,bb/ m b 

Figure 3: BR of B — ► 7TT + T~as a function of £jy bb/ m b f° r £jv rT = 10 m T and |r t b| < 1. Here Bi? is 
restricted in the region between solid (dashed) curves for Cf f > {Cj Sf < 0). Small dashed straight line 
represents the SM prediction. 




Figure 4: The same as Fig.(^j), but for r t b > 1 with £jy TT = m T . 



are extremely large. These figures explicitly show the size of the NHB effects on the exclusive 
B — ► ttt + t~ decay. 

From now on in figures we plot, the regions bounded by the solid (red) curves represent the 
> case and the regions bounded by the dashed (green) curves represent the Cy < case 
while the small dashed (blue) curves are for the SM predictions for the relevant observable. 





BR(B -» ttt+t^) 




\V cb \ = 0.037 \V cb \ = 0.043 


(0.3; 0.34) 
(-0.3; 0.34) 
(-0.07; 0.34) 


0.69 x 10" 8 0.93 x 10" 8 
0.60 x lO" 8 0.81 x 10~ 8 
0.62 x 10" 8 0.83 x 10" 8 



Table 3: The values of the total branching ratio for B — > ttt + t decay in the SM, at three different sets of 
the Wolfenstein parameters (p; 77). 
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Figure 5: BR of B — ► 7rr + T~as a function of £jy TT for ^ bb = 40 mj, and \r t b\ < 1. Here £>i? is restricted 
in the region between solid (dashed) curves for > (Cy < 0). Small dashed straight line represents 
the SM prediction. 
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Figure 6: Differential Afb for £> — * 7rr + r as a function of s for £jy bb = 40 m& and ^ TT = 10 m T , in 
case of the ratio \r t b\ < 1. Here differential Afb is restricted in the region between solid (dashed) curves 
for C e 7 ff > (Cf f < 0). 

We present the dependence of the BR on the parameter fjy bb / m b m Fig s -(§"0)' where the first 
one is for \r t b\ < 1 and the latter for r#, > 1. Our prediction for the BR of the B — > 7rr + r _ decay 
in the SM including the long distance effects is 

BR(B -> vrr+r") = 0.76 x 10 -8 . (3.7) 

We also give the SM values of the total branching ratio for B — > 7rr + T _ decay at three different 
sets of the Wolfenstein parameters (p;r/) in table |[ As seen from Fig. (|3]> , for \rn\ < 1, the 
Cy < case almost coincides with the SM prediction. However, for C7 > 0, we observe an 
enhancement which is 1.5 — 2 times of the SM prediction; but this enhancement decreases with 
the increasing values of the Cjv6b/ m & parameter. In case of r t b > 1 (Fig. (|j)), extremely large 
enhancement, 2 — 3 orders larger compared the SM case, is reached for Cj < case. 

We plot the dependence of the BR on the parameter £jy Tr in Fig. (|J) for \r t b\ < 1 . From this 
figure, we again observe an enhancement as in the fjy bb / m b dependence and this is the contribution 
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Figure 7: The same as Fig.(g), but for r tb > 1 with bb — 0.1 m b and TT — m T . 
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Figure 8: Afb for B — > ttt + t as a function of £5 rr for £jy ftb = 40 and |r t fc| < 1. Here ^4fb is 
restricted in the region between solid (dashed) curves for Cf f > (Cf f < 0). 
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$,rr (GeV) 

Figure 9: The same as Fig.(^j), but for r tb > 1 with |jy bb = 0.1 m b . 

due to the NHB effects. However, the behavior of this dependence is opposite to that of |^ bb /m b 
dependence: the BR increases with the increasing values of the ^ TT . The SM prediction again 
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lies in the region bounded by the Cy < case. 

The dependence of the differential Afb on the invariant dilepton mass s for the B — > -kt + t~ decay 
is presented in Fig. (Q) (Fig. ([7])) for |r#,| < 1 (r^ > 1) case. Since Afb arises in the 2HDM 
only when the NHB effects are taken into account, it provides a good probe to test these effects. 
For \r t f,\ < 1, although Afb{s) is very small for Cj < case, it is considerably enhanced for 
Cy > case. For r t b > 1, Afb(s) for < and C7 > cases completely coincide and 
its magnitude is one order smaller than the |j4f\b(s)| for \r t b\ < 1 case. 

Figs.(|J) and @ are devoted to the £jy TT dependence of the Afb of the B — * 7rr + r _ decay for 
\r t b\ < 1 and r t b > 1 cases, respectively. As can be observed from Fig.(||), Afb is quite sensitive 
to the parameter £jy Tr especially for Cy > 0. It can reach 10% for £% TT ~ 60. For r t6 > 1 
case shown in Fig.(^), the Cy * < and C^* > cases completely coincide and l-Ap^l decreases 
with the increasing values of the fjy TT . In addition, its value is one order smaller than the |r#,| < 1 
case. 



3.2 Numerical results of the exclusive B — ► decay 

In our numerical calculation for .B — » /9£ + ^~decay, we use three parameter fit of the light-cone 
QCD sum rule [ 30 ] which can be written in the following form 

F (l 2 ) = ~\ 2 / 2~TZ~T~2~/ — 2T2 

1 — clf q /m B + OF{q /m B ) A 

where the values of the parameters F(0), clf and bp are given in table @. The form factors Aq 
and ^3 can be found from the following parametrization, 

4, = A, - 



m p m b 



A, = H^R Al _ !!^RA 2 . (3.9) 




Figure 10: dBR/ds for B — > pr+r as a function of s for ^ bb = 40 mf, and £jy rr = 10 m r , in case of 
the ratio \r t b\ < 1. Here dBR/ds is restricted in the region between solid (dashed) curves for C e 7 fS > 
(C7 < 0. Small dashed curve represents the SM prediction. 
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Figure 11: The same as Fig.(fTo|), but for r t b > 1 with £jy bb = 0.1 nit, and £jy TT = m T . 
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Figure 12: BR for B — > pr+r^as a function of £jy bb /mb for £^ rT = 10 m r and \rtb\ < 1- Here £>i? is 
restricted in the region between solid (dashed) curves for Cy ff > (Cj ff < 0). Small dashed straight line 
represents the SM prediction. 



We first consider the dependence of dBR/ds on the invariant dilepton mass s for the B — > 
pr + r~decay. This is plotted in Fig.(pTj|) for CjvM) = ^0 m b an ^ Ijvtt = ^ mr ' * n case °^ tne 
ratio | | < 1 by taking into account the long distance effects. We conclude from this graph 
that, dBR/ds almost coincides with the SM result for Cj < case, while for > it is 
considerably enhanced. As for the r t b > 1 case (Fig. (|Tl|))where we take £® bb = 0.1 m& and 
£n tt = m T , we observe an enhancement of one order as compared with the \r t b\ < 1 and also the 
SM cases. Here Cf f > and C e 7 Sf < cases completely coincide. 

The dependence of the BR on one of the free parameters of the model III, |^ bb /m b , is given 
in Figs. ( |l2| ) and (13) for \rtb\ < 1 and r t b > 1, respectively. Our prediction for the BR of the 
B — ► /9T + r~decay in the SM including the long distance effects is 

BR(B -> pr+r") = 0.94 x 10~ 8 . (3.10) 

We also give the SM values of the total branching ratio together with Afb values for B — > 
pT + T~ decay at three different sets of the Wolfenstein parameters (p; rj) in table f|. As seen from 
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Figure 13: The same as Fig.([l2|), but for r t b > 1 with £jy TT = m T . 
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Figure 14: BR for B — ► pr+r^as a function of |^ TT for bb — 40 m b and |rt&| < 1. Here Si? is 
restricted in the region between solid (dashed) curves for Cf f > (Cj ff < 0). Small dashed straight line 
represents the SM prediction. 





BR(B -> pr+r") 


A FB {B -► pr+T-) 


{p; v) 


\V cb \ = 0.037 


\Vcb\ = 0.043 




(0.3; 0.34) 


0.94 x 10" 8 


1.26 x 10" 8 


-0.17 


(-0.3; 0.34) 


0.72 x 10~ 8 


0.97 x 10" 8 


-0.20 


(-0.07; 0.34) 


0.76 x 10~ 8 


1.02 x 10~ 8 


-0.19 



Table 4: The values of the total branching ratio and the forward-backward asymmetry for B — > pr+r decay 
in the SM, at three different sets of the Wolfenstein parameters (p; rf). 

Fig. (|l2|), for \rtb\ < 1 where we take ^ TT = 10 m T , the Cy < case coincides with the SM 
prediction. When > 0, however, the BR is enhanced by 2.5 — 5 times of the SM prediction; 
but this enhancement decreases with the increasing values of the fjy bb / m b parameter. For r tb > 1, 
we take fjy rr = m T and observe an enhancement for both Cy < and > cases. For 
C 7 e// < 0, it is one order larger than the SM value while for Cj^ > 0, the order of enhancement 
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Figure 15: Afb(s) for B — ► pr + r as a function of s for £® b6 = 40 m& and £^ rT = 10 m T , in case of 
the ratio \r t b\ < 1. Here A^b(s) is restricted in the region between solid (dashed) curves for Cj ff > 
(Cy < 0). Small dashed curve represents the SM prediction. 
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Figure 16: The same as Fig.(p3[), but for r tb > 1 with £®.bb = 0.1 m& and £jy TT 



is the same as that in |r#,| < 1 case. 

We plot the dependence of the BR on £jy , the other free parameter of model III, in Fig. 
(14) for |r ib | < 1. Here, we take £5 bb = 40m ft and see that the BR is not sensitive to TT for 



N,bb 

< and it is almost the same as the the SM prediction. However for > 0, the BR is 
quite sensitive to fjy TT and increases as it increases. 

The dependence of Afb(s) of the B — > pr + T~decay on the invariant dilepton mass s is 
represented in Fig.(|l~5|) (Fig.([Tr^)) for \r tb \ < 1 (r#, > 1) case. For \r tb \ < 1, there is an enhance- 
ment on |Af#(s)| for > 0, while for C% < it is almost the same as the SM prediction. 
For r t b > 1, all Model III predictions for Afb(s) almost coincide but with a flip in the sign as 
compared to the SM prediction. 

Finally we present the dependence of the Aps on the ^Ntt parameter. Our prediction for the 
Afb of the B — > pr + r~decay in the SM is 

A FB (B -> P t + t-) = -0.193 (3.11) 

(See also table f|.) In Fig. (|l7|), this dependence is plotted for the ratio \r t b\ < 1 and |^ bb = 40 m b . 
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Figure 18: The same as Fig.(|l7|), but for r t b > 1 with £jy bb = 0.1 mj. 



Although the C% ff < case almost coincides with the SM value, there is an enhancement for 

cD . 

*N,TT 



\Afb \ up to the 50% of the SM value for the moderate values of the £5 rr parameter, for > 



case. For r t b > 1 case where we take £® bb = 0.1 my,, \ Aps\ can reach at most half of the SM value 
and drops to zero for large values of £jy TT - Here, Cy < and > cases almost coincide. 

Finally, we would like to comment briefly about the NHB effects on the CP violating asymme- 
try, Acp, for B — ► 7r£ + ^~and B — >■ p£ + £~decays. As pointed out before @]-[@], in the SM there is 
a considerable Acp in the partial rates for these decays because all three CKM factors contributing 
are at the same order. Further, the 2HDM contributions to Acp have been investigated in [Q, 10] 



and it is shown that since charged Higgs contributions give rise to constructive interference to the 
SM result, Acp decreases while the BR increases for B — ► 7r^ + £~and B — ► p^+^decays in the 
2HDM. We expect that including the NHB effects will further decrease magnitude of Acp with 
respect to its value without NHB effects. To see this, consider Aqp between B — ► M£ + £~ and 
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B — » M£ + £ decays for M = it, p, which can be written as [20, 31] 



A 



CP 



f dsIm{C e 7 ff )Im{Cl ff )F 
/ ds A 



where T is a function of various form factors for the decays we consider and A is proportional to 



the one given by Eq.(2.16) and (2.26) for ir and p, respectively. As can be seen from the equation 
above, the numerator of the Acp ratio is free from the NHB contributions while the denominator 
gets this additional contribution so that magnitude of Acp will decrease with the inclusion of the 
NHB contributions. 



3.3 Conclusion 

In this paper we have investigated the physical observables, BR and Afb, related to the exclusive 
B — ► ir£ + £~a.nd B — ► p£ + £~decays in the general 2HDM including the NHB effects. We have 
found that NHB effects are quite sizable, leading to considerable enhancements on these physical 
observables. An experimental observation of the Afb in the B — > 7r£ + £~ decay, which is absent 
in the SM, would be a very powerful and direct test of the 2HDM and the existence of NHB. In 
conclusion we say that the exclusive B — ► ir£ + £~a.nd B — > p£ + ^ _ decays provide very useful 
testing ground for the new physics beyond the SM. 



A. The operator basis 



The operator basis in the 2HDM for our process is [15, B2L 33] 



01 = (dLalfiCL^icL^bLa), 

2 = (d~La'ynCLa){cLf3'y fJ 'b L f3), 

3 = (d La i^b La ) Y (qL^qip), 

q=u,d,s,c,b 

04 = (d~L a Jf,b Lf 3) Y (QLpl^qLa), 

q=u,d,s,c,b 

q=u,d,s,c,b 

6 = (dLal^b L/3 ) Y (qRP^QRa), 

q=u,d,s,c,b 

g _ 

07 = Y^2d a a l _ lu (m b R + m s L)b a J : ' 11 ', 



j^d a T^a^ u (m b R + m s L)bpQ 
(d La ^b La )(£^£) , 



01 = (d La -i^u L f3)(u LI 3^b La ) 

02 = {dL a l^u La )(u L ^b L(i ) 



afiu 



16vr 2 

e 
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Q2 = ^(dlb%)(t l5 £), 

q=u,d,s,c,b 
q=u,d,s,c,b 

q=u,d,s,c,b 

q=u,d,s,c,b 
q=u,d,s,c,b 

Q9 = ^(dl^b P R ) Y, (fl£«V«S). 

010 = ^(^^65) E (flS^flZ), (A.i) 

where a and /3 are 5t7(3) colour indices and and Q^ v are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 



B. The initial values of the Wilson coefficients. 

The initial values of the Wilson coefficients for the relevant process in the SM are [ ]32| ] 

-rSM 



C* M (m w ) = 1 , 

C 7 SM (m w ) = ]f - 2X l \nx t + " 8X * " 5X * + 7Xt 



4(x t - l) 4 11 2A(x t - 1 



,3 



3xf — xf + 5xf + 2x t 



Ci M (m w ) = - A ,f Xt ^ d lnx t + _ t]i 

sin A &w sin Ow 9 



4(xt - l) 4 8(x t - 1) 

\ 1 \ 1 — 4sin 2 6w 4 
Ci , (m w ) = r^— S(s t ) + — -C(s t ) - ^fo) + 



cf A W) = -in^) - c i x t)) , 

sm fvi/ 

Cgf(m^) = * = !,.., 10 (B.l) 
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and for the additional part due to charged Higgs bosons are 



Ci%(™ W ) = 0, 
Cf (m w ) = Y 2 F 1 (y t ) + XY F 2 (y t ) , 
C^(m w ) = Y 2 Gx(y t ) + XY G 2 (y t ) , 
C?(m w ) = Y 2 H 1 (y t ) , 

C? (m w ) = Y 2 Ltiyt) , (B.2) 



where 



Y = — (&, tt + tN,tM) ■ (B-3) 

Note that the results for model I and II can be obtained from model III by the following substitu- 
tions: 

Y -> cot /3 , XY - cot 2 p for model I 

Y -> cot /? , 17^1 /or mode/ /J . 



The NHB effects bring new operators and the corresponding Wilson coefficients read as [11] 

CN,rT^N,tt) 3m byt{Q5(yt)zA - &i(z A ,y t )) 



CQ 2 ((£,N,tt 

Cq 2 ((£n,u)' 2 ) 



32Tr 2 m 2 AO m t Qi{zA,yt)&5(yt) 



A°/(tU \2\ _ ^N,TT\^N,tt) $N,bb 



32TT 2 m 2 A0 \®i(zA,yt)®i{zA,yt)&5(yt) 

(y t (@i(z A , y t ) - Q 5 (y t )(xy + z A )) - 2Q 1 (z A , Vt)*s(Vt) H ^^k ) 

<vi{ZA,yt) 



r A°,pU \ _ g2 ^N,TT^N,tt m b x t ( 2 xyx t + 2z A , r z A Q 5 (x t ) 1 

Q 2 ^,ftJ- 647r 2 m 2 omt y @5 (x t ) e i (z A ,X t ) ^Gtiz^Xt) 1 

- xyx t y t { — - l ^ Xt ^ Vt / ZA ~ 1 ) ~ ( 1 + x )yt) 

(0 6 - (x - y)(x t - y t ))(@3{z A ) + (x - y){x t - yt)z A ) 

x(y t + x t {l - yt/zA)) ~2y t X 



e 6 e 3 (z A ) 



-) 



■At LI CU 

G Q^N,bb) ~ 64vr 2 m 2 



9 Z CN,TT^N,bb A _ x 2 t y t + 2y(x - l)x t y t - z A {x 2 + e 
64ir 2 m 2 A0 V e 3 (z A ) 



+ x 2 (l-y t /z A ) | 2 z A 6 



©6 Q 2 (za,x) 
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C Ql MN,U) ) ~ 647r 2 m 2 



9 2 {iN,tt? m bm T fx t (l - 2y)y t (-1 + 2 cos 2 6 W )(-1 + x + y)y t 



2 

H o m t 



e 5 ( yt ) 



+ 



cos 2 w e 4 (yt) 



zg(Qi(zg, j/t)a:yt + cos 2 W (-2x 2 (-l + x t )yy t 2 + xx t yy 2 - 6 8 z g )) \ 



cos 2 8 w @ 1 (z H ,y t )@ 7 



riH ( fu \ _ 9 2 ^N^N,bb m r ( (-1 + 2 cos 2 9 w )y t _ xtyt Xtytjxy - z H ) 
q 1 ^n,u)- 647r 2 m | omt I cos 2^e 4(yt) ' e 5 (y t ) @i(z H ,yt) 



(-1 + 2 cos 2 6 w )y t z H 
+ 2~TT\ 2xt ln 

COS z dw^>7 



&i(zH,yt) 



(B.4) 



,//' i g 4 m b m T x t ( ^ | (-l + 2x)x f 

' ff0 mt I ©5(^t) + y(i - x t ) 



CQl (ff } " ~128vr 2 m 2 



+ 



2x t (-l + (2 + x f )y) 

e 5 (x t ) 



4 cos 2 gwj-1 + x + y) + x t (x + y) x t (x(x t (y - 2z H ) - ^z H ) + 2z H ) 

cos 2 9 w @ 4 (x t ) &l(ZH,X t ) 
y t ((-l + x)x t z H + cos 2 <V((3x - y)z H + x t (2y(x - 1) - z#(2 - 3x - y)))) 



cos 2 0w(@3(zh) + » (art - yt)z#) 



+2 (x* In 



©i(-Zif,ar t ) 



+ ln 



x(y t - x t )z H - @ 3 (z H ) 
(G 5 (x t ) +y(l - x t )y t z H 



C h Q \{(f N , tt ?) = - 



&rM,tt) 3 ™ b yt 



ei(z h ,y t ){2y - 1) + e 5 (y t )(2x - l)z h ) 



327r 2 m 2 mtei(z/ l ,y t )e 5 (yt) 

^N,TT^N,tt) 2 ( (®b(yt)zh(yt - l)(ar + y- l) - ©i(^,yt)(©5(yt) +yt) 



32vr 2 m 2 o 



©i(^)©5 (yt) 



21n 



Zh®5(yt) 
Qi(zh) 



9 2 ^ T JN,tt m bXt ( 2(-l + (2 + x t )y) x t (x - l)(y t - z h ) 



+ 



647r 2 m 2 mi ^ ©5 (art) 
x(x t (y - 2z ft ) - 4z fe ) + 2^ 



© 2 (^) 
(1 + x)y t z h 



+ 21n 



@l(z h ,x t ) 



@i{z h ,x t ) 



xyx t y t + z h ((x - y)(x t - y t ) - ©e) 



©9 + yt^fe((ar - y)(x t - yt) - © 6 )(2x - 1) x(y t z fe + x t (z ft - y t )) - 2y t z h 



^@6(©6 - (ar - y){x t - y t )) 



© 2 (^) 



/ff \ 9 ^N,rr^N,bb ( yx t y t (xx 2 (y t - z h ) + @ 6 z h (x - 2)) 
°Qil?Ar,fefej r-r^ — o — I 77-^ — — rzm 



64vr 2 



m 



h° 



&6_ 

xtyt 



+ 2 In 



x t ytzh 

@2(z h ) 
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where 



A) = —(—1 + y — yX)uj — x(y\ + lu — uj\) 
9 2 (w) = (x t + y(l - x t ))y t w - xx t {yy t + (y t - 1)cj) 
@' 2 (lo) = Q 2 (u,x t <-» y t ) 

@ 3 (u) = + y) - y)y t w + xx t (yy t + w(-l + y t )) 

04 (cj) = 1 — X + XCJ 

e 5 (A) = x + A(l-x) 

©6 = (xt + 2/(1 - £f))yt + xx t (l - 2/t) 

@7 = (y(yt - 1) - ytW + x(yy t + (2/t - 1)2//) (B.5) 
9 8 = yt(2x 2 (l + z t )(jft - 1) + a?t(y(l - 2/t) + yt) + x(2(l - y + yt) 

+ x t (l-2y(l-y t )-3y t ))) 
©9 = -^t (-1 + x + y)(-yt + sc(2yt - l))(yt - z h ) - x t y t z h (x(l + 2x) - 2y) 

+ y\( x t(x 2 - y(l - »)) + (1 + x)(x - y)^) 



and 



act 



771 



yt 



ir 



777 



2ft 



mi 



mi 



H° 



777 



h° 



777 



The explicit forms of the functions F l(2 )(yt), G l{2 )(y t ), Hi(y t ) and Li(y 4 ) in Eq.( 
as 

yf) , y t 2 (3y t -2) 



Fi{y t ) 
F 2 (y t ) 
Gi(y t ) 
G 2 (y t ) 
Hi(y t ) 

Li(vt) 



+ 



yt(7 - 5y t - 

72(^-1)3 ' 12(^-1)4 
yt(%t - 3) yt(-3y 4 + 2) 

Inyt 



In y* 



+ 



12(y t - l) 2 6(y t - 1) 

yt(-y 2 + % + 2) 



24(y 4 - 1)= 



+ 



"2/t 



4(y t - 1) 4 



In yt , 



+ 



yt 



ytjyt - 3) 

4(y t -l)2 1 2(y t -l)3 
1 — Asin 2 6\Y xyt 



In yt , 



sin 2 8 



yt 



y t - 1 (y t - l) 2 
47y 2 - 79y t + 38 3y 3 - % + 4 
18(yt - l) 4 
1 



In y t 



W8(y t - If 



In y t 



1 



xyt 



sin 2 0w 8 



+ 



yt - i (yt - i) : 



In y t 



are given 



(B.6) 



Finally, the initial values of the coefficients in the model III are 



C 



2HDM 



(m w ) = c! M {m w ) + C?(mw), 



C, 



2HDM 



(m w ) 



dx 



l-x 



dy (<((&) 2 ) + <(&) + <(y 4 ) + <((e5U) s ) 
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J o ■/ o 

(~i2HDM j \ m b ln 2HDM lrr ^ \, n 2HDMi^ \\ 

rn^sm Uw 

L Qi \ m W) = . 2a y L Qi \ m w)-^Q 2 ymw)) 

rri£ sm u\y 

C%* DM (m w ) = , i = 5,...,10. (B.7) 

Here, we present Cq 1 and Cq 2 in terms of the Feynman parameters x and y since the integrated 
results are extremely large. Using these initial values, we can calculate the coefficients Cf HDM (/i) 
and CqJ dm (fi) at any lower scale in the effective theory with five quarks, namely u,c,d,s,b 
similar to the SM case Q-Q. 

The Wilson coefficients playing the essential role in this process are Cj HDM (//), Cg HDM (//), 
CIq DM (pl), CqI dm (pl) and Cq^ dm (//). For completeness, in the following we give their explicit 
expressions. 

c?"(m) = cr DM {n) + Qd (ct HDM (n) + N e c! HDM (ri) , 

where the LO QCD corrected Wilson coefficient Cj° ,2HDM (fi) is given by 

c LO,2HDM^ = ^S^HDM^ + (g/3) ^14/23 _ ^16/23 j ^HDM ^ 

8 

+ Ci HDU (m w )J2hiV ai , (B.8) 

i=l 

and 7] = a s {mw) / ch s {ijl), hi and ctj are the numbers which appear during the evaluation [|36|]. 

Cg (fi) contains a perturbative part and a part coming from LD effects due to conversion of 
the real cc into lepton pair £ + l~ : 

C e 9 ff (fi) = C p 9 ert (n) + Y reson (s) , (B.9) 

where 

+ h{z, s)[3Ci(/i) + C 2 (/i) + 3C 3 (^) + C 4 (m) + 3C 5 (^) + C 6 (/x) 
+ A u (3d + C 2 )] - ^(1, s) (4C 3 f» + 4C 4 (^) + 3C 5 (fi) + C 6 (jjl)) 

- ~h(0, s) [C 3 (m) + 3C 4 (m) - A u (6Ci( M ) + 2C7 2 (//))] (B.10) 

+ ~ (3C 3 (/u) + C 4 (m) + 3C 5 Gu) + C 6 (m)) , 



and 



Y (?) - — ——— k V nm - e+r ^ 
"em 9 2 - m Vi + im Vi T Vi 

x [(3Ci(//) + C 2 (M) + 3C 3 (^) + C 4 ( M ) + 3C 6 (aO + C 6 (m)) 

+ A u (3Ci( M ) + C 2 (//))]. (B.ll) 
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In Eq.( p.9| ), the functions h(u, s) are given by 



h(u, s) 



■— in — 
9 n 



- Inu H — - H — x 
9 27 9 



9 



(2 + s)|l - x 




27T ) , for x 



Au 2 



< 1 



for x = ^- > 1, 



ft U, s = In In s H — Z7r , 

V ' 7 27 9 /x 9 9 ' 



(B.12) 



(B.13) 



with u = S* 6 . The phenomenological parameter k in Eq. ( [B.llD is taken as 2.3. In Eqs. (B.10) and 
(B.ll), the contributions of the coefficients C±(fi), Ce(fi) are due to the operator mixing. 



Finally, the Wilson coefficients Cq^h) and Cq 2 (h) are given by [15] 

C Q ^)=V~ 12/23 C Ql (m w ) : i = l,2. 



(B.14) 
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